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Repeated studies of the cosmic microwave background (CMB) based on data from the Wilkinson 
Microwave Anisotropy Probe have revealed an apparent asymmetry in the distribution of tempera- 
ture fluctuations over the celestial sphere. The studies indicate that the amplitudes of temperature 
fluctuations are higher in one hemisphere than in the other. We consider whether this asymmetry 
could originate from a large scale inhomogeneity in the gravitational field enclosing the present 
Hubble volume. We examine what effect the presence of an inhomogeneity in the gravitational field 
of size larger than the present Hubble radius would have on the temperature distribution of the CMB 
and start eliciting its observational signature in the CMB power spectrum. The covariance function 
(a* m a;' m /) contains, in addition to the diagonal Ci entries of the conventional CMB temperature 
anisostropy power spectrum, non-diagonal entries. We find that specific non-diagonal entries of the 
covariance function are sensitive to the strength of the inhomogeneity, while the diagonal C'i entries 
are not. These non-diagonal entries, which are not present in the case of a homogeneous background 
geometry, are observational signatures of a large-scale inhomogeneity in the background geometry 
of the universe. Furthermore, we find that an inhomogeneity in the gravitational potential of super- 
Hubble size would yield a power asymmetry in the CMB with maximal asymmetry at an angle of 90° 
, to the CMB dipole axis. The axis of the CMB power asymmetry was recently estimated by Eriksen 

et. al. to be located at angles between 83° and 96° to the CMB dipole axis, which is consistent with 
the prediction of our model. This implies that the location of the observed power asymmetry in the 
^ C"| CMB sky could be accounted for by a large-scale inhomogeneity in the gravitational field enclosing 

O !■ the present Hubble volume. 
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PACS numbers: 98.80.-k, 98.80.Es, 98.70.Vc 



^ . I. INTRODUCTION 

> . 

Precision measurement and analysis of the cosmic microwave background (CMB) over the past two decades has 
exposed an amazingly detailed picture of the early universe, a picture in which not only the processes and structure 
of the early universe is imprinted, but also the signature of any subsequent gravitational processing of this structure. 
This makes the CMB an extraordinarily rich and powerful tool for probing and depicting physical processes in the 
early universe. 

Hence, analysis of the CMB has uncovered evidence of how the observed structure in the universe has emerged: In the 
primordial universe, stochastic processes in the form of quantum fluctuations in primordial fields drove fluctuations 
• • i in the gravitational field which subsequently gave rise to fluctuations in the matter and photon densities. These 
primordial density fluctuations are today imprinted in the temperature fluctuations in the CMB, and the properties of 
the primordial processes are encoded in the statistical properties of the CMB. In order to understand the primordial 
universe, it is therefore of utmost importance to understand the statistical properties of the CMB and decode its 
statistical information. 

Full-sky analysis of data from the Wilkinson Microwave Anisotropy Probe (WMAP) are in good agreement with 
the predictions of simple inflationary models that prescribe a nearly scale free, Gaussian primordial power spectrum 
and a universe that is flat and homogeneous at large scales 0, HJ . The CMB spectrum expected in such a universe 
is statistically isotropic, implying that the expectation value of CMB observables are rotationally invariant. There 
is, however, an apparent asymmetry in the power of CMB flucutations between two different hemispheres J3041| 
of the CMB sky. Various explanations of the anomaly have been proposed. New inflationary models [7-9] may 
yield asymmetric power spectra, which may account for the CMB power asymmetry. Other authors have invoked 
conventional physics, such as non-linear gravitational effects of local inhomogeneities [10] , or dust-filled voids in the 
local universe 
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Here, we are taking a simplistic phenomenological perspective where we consider whether the presence of a large 
scale inhomogeneity in the gravitational field enclosing the present Hubble volume could account for the power 
asymmetry|23|. At this stage, we are not concerned with how the inhomogeneity came about, but it is natural to 
assume that it has primordial origin. However, independent of its origin, a gravitational inhomogeneity in the form 
of a large scale perturbation to the gravitational potential could possibly manifest itself through a visible asymmetry 
in the CMB power spectrum. If successful, we may be able to capture the essentials of the physical phenomena under 
study in a rather model independent way, without excluding relevant models that may account for the phenomena. 

If we were to detect large-scale asymmetries in the CMB temperature distribution, the determination of their origin, 
whether they originate from asymmetries in the primordial power spectrum or from large scale inhomogeneities in 
the gravitational field, would in any case provide us with valuable additional insight into the nature of the processes 
in the primordial universe. 

Given an early inflationary universe, large scale fluctuations in the gravitational field will directly depict primordial 
processes, because as they leave the horizon early in the inflationary era, large scale scalar perturbations to the 
gravitational field freeze. Moreover, they would stay more or less constant until now, only decaying slightly throughout 
the history of the universe [lH • 

How could such large-scale perturbations to the gravitational field - of size larger than the observable universe - 
be observed? This is the first question that we would like to answer in this paper. Under the assumption that our 
observable universe is immersed in a scalar metric perturbation of super-Hubble size, the aim of this paper is to 
determine what effect the presence of such a perturbation would have on the CMB temperature and start to elicit its 
observational signature in the CMB spectrum. 

Given the puzzling evidence of a hemispherical power asymmetry in the CMB a natrual follow-up question is: 
Could an hemispherical power asymmetry in the CMB be caused by a large-scale inhomogeneity in the gravitational 
field of super-Hubble size? Posing these questions, we make no assumptions about the origin of his perturbation 
or any physical mechanism that might have caused it. Our goal is simply to determine what effect the presence of 
such a perturbation to the gravitational potential would have on the observed CMB temperature and to elicit its 
observational signature in the CMB spectrum. 

Naively, one could perhaps attempt to answer these questions by looking for an axis of maximal anisostropy in the 
CMB spectrum [3] and correlate the apparent CMB anisotropy with the expected anisotropy from a super-Hubble 
perturbation to the gravitational potential. There is, however, already such an axis of maximal anisotropy: the CMB 
dipole, which originates from our motion relative to the CMB rest frame and by far dominates the CMB anisotropy 
spectrum. We will therefore have to look for the signature of a secondary anisotropy axis in the residual temperature 
map after first removing the CMB dipole. 

This paper starts by computing the expected temperature anisotropy from a super-Hubble perturbation to the 
gravitational potential, as seen from the CMB rest frame. Then, we transform this result to a moving frame of 
observation, before we factor out the best fit dipole. Treating the resulting residual as a lowest order perturbation 
to the temperature distribution expected from a standard, homogeneous background cosmology, we can therefore 
factor the residual into the temperature distribution of an arbitrary homogeneous background cosmology. By doing 
a standard multipole expansion of the temperature distribution, we are finally able to obtain the signature effect of 
the super-Hubble perturbation on the CMB spectrum. 

II. SUPER-HUBBLE PERTURBATIONS 

Let us briefly review the large scale solutions to the first order perturbation equations for a ACDM cosmology. 
More details can be found in Appendix [AJ and in (l2| . In conformal Newtonian gauge, the space-time metric takes 
the form 

ds 2 = a 2 (f]) (- (I + 2^(77, x)) Arj 1 + (f + 2$(?7, x)) <%dxMx J ) , (1) 

assuming implicit summation of spatial indices, rj is the conformal time, while ^ and $ are the two scalar metric 
perturbations. In this section, we will restrict our attention to the evolution of the gravitational potentials $ and ty, 
and will therefore discard any variable without significant contribution to their evolution. 

Before we continue, let us define what we mean by a super-Hubble perturbation. As we primarily will work in the 
spatial domain, we need a definition of super-Hubble perturbations in the spatial domain that is compatible with the 
conventional definition in the frequency domain. In the frequency domain, a super-Hubble perturbation is a mode 
with wave length much larger than the comoving Hubble radius W , where TL = a/a. A super-Hubble mode with 
wave vector ki will therefore have k 2 << W 2 . Similarly, for any perturbation variable X, we define a spatial domain 
super-Hubble perturbation to be any spatial perturbation with spatial derivatives that satisfy |9iA/A| << % and 
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\d 2 X/ X\ « H 2 . By applying this definition to a Fourier mode ~ e tk ' x , we see that the definitions of spatial-domain 
and frequency domain super-Hubble perturbations are consistent. 

At super-Hubble scales, the temporal dependence of the perturbation variables can be separated from the spatial 
dependence, allowing us to write the potential as a product of a fixed spatial perturbation and a time-dependent 
factor: 

$ ~ ~ $0?)$o(» (2) 

where <f>o(x) is the primordial spatial perturbation. Now, since $ is a super-Hubble perturbation, we have 
W~ 1 di$>o << $o- Consequently, within the present Hubble volume, higher-order derivatives of $o can be ignored, 
and $o may be approximated by its first order Taylor expansion around the position x of the observer: 

$ (a;) ~ <Z> (x Q ) + di<Z> (x ) (x* - xi) (3) 

The spatial coordinates {x 1 } are arbitrary, and we need to fix them in order to proceed. Let us choose cartesian 
spatial coordinates x, y, z that make the potential, and therefore the metric of eq. ([T]) manifestly axially symmetric 
to first order in di$o(xo). Align the z axis with the gradien V^o^o) an d let the origin be the observer's position. 
Then the linear expansion of the potential in eq. ([3]) takes the form 

$ (z)~-fc(z-z ) (4) 

1/2 

where the constant k = (5 13 di<&o(xo)dj$o{xo)) and zq is a constant. With this approximation, the metric of eq. 
(fTj) now takes an axially symmetric form: 

ds 2 = a 2 (r]) (- (1 + 2^(77, z)) Arj 2 + (1 + 2$(r?, z)) (dx 2 + dy 2 + dz 2 )) (5) 

In the following, we will use the space-time metric of eq. ([5]) as the background metric. This background can 
be treated as a perturbation to the homogeneous background ds 2 = a 2 {ri)(—dri 2 + 8^dx l dx^). To first order in the 
perturbation variables, this pertubation evolves independently of other perturbations, a result of the perturbation 
equations being linear. We are therefore able to treat the effect of a super-Hubble metric perturbation independently 
of other effects. 



III. THE EFFECT OF A SUPER-HUBBLE METRIC PERTURBATION ON THE ANGULAR CMB 

TEMPERATURE DISTRIBUTION 

In order to compute the CMB temperature perturbation arising from a super-Hubble metric perturbation, we need 
to solve the equations of motion for photons arriving from an arbitrary direction (0,ip). In the following, we will 
assume a space-time metric on the form 

ds 2 = a 2 (rj) (-A{t], z)d V 2 + A(r], z)^ 1 (dx 2 + dy 2 + dz 2 )) , (6) 

where A(r), z) = 1 + 2\P(?7, z), which is an idealization of the general case, but reduces to eq. §5§ to first order in 3> 
in the case "J = — $. 

Once we have computed the trajectories of free-streaming photons arriving from an arbitrary direction, we will 
be able to compute the gravitational potential at the time of emission. Assuming CMB photons were emitted 
simultaneously at a uniform temperature at recombination, the current CMB temperature will primarily depend on 
the gravitational potential at the time and place of emission as well as any subsequent development of the gravitational 
potential along the paths of the CMB photons. 

A. Computing null geodesies 

The first item on our agenda is to compute the trajectories of photons moving in the inhomogeneous background 
defined by eq. These trajectories are null geodesies, determined by the geodesic equations. The geodesic equations 
are derived by minimizing the action S = J dX-^g^^^-^- for a massless particle moving in the background of eq. 
©. Null geodesies, as well as certain geometric quantities, like angles, are invariant under conformal transformations. 
We therefore choose to compute the null geodesies in a different conformal frame, using the conformally related metric 



ds 2 = -A{ri, z)dr] 2 + A(ri, z)' 1 (dx 2 + dy 2 + dz 2 ) 



(7) 
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The reason for this is that the equations of motion simplify considerably in this metric. Notice that, even if the 
null geodesies and certain other quantities are invariant under conformal transformations, most physical properties of 
the photons depend on the conformal frame. We will therefore take care to use the physical spacetime metric of eq. 
([5]) when computing conformally variant physical properties of the photons, such as the photon energy. 

Two constants of motion, v x and v y , arise because coordinates x and y are cyclic coordinates of the action. The 
equations of motion for coordinates x and y are, using overdots to denote differentiation with respect to the affinc 
parameter A: 

A' 1 x = v x ,A- 1 y^v y (8) 
Furthermore, the equation for the time coordinate rj is 

Finally, the equation for the z coordinate becomes 

1 OA 8 A 2 

z- - — zr, + A—v T 2 = 10 
A or\ oz 

where the constant vt is defined by vt 2 = v x 2 + v y 2 . These equations of motion are generally not solvable by 
analytical means. However, we will introduce an idealized model that allows analytical treatment. The model is 
motivated by the observation that, as shown in Appendix large scale perturbations to the gravitational potential 
stay more or less constant until the universe ceases to be matter dominated. 

The idealized model that will be used from now on asserts a time- independent gravitational potential with = 
const = 1 in the metric of eq. ([5]). In this case, the metric is conformally static, with 

A(z) = 1 + 2k(z - z ) = A + 2kz (11) 

Notice that this is an unrealistic approximation at all but the largest scales. At large, super-Hubble scales, still 
being a very crude approximation, it is nevertheless a useful idealization that grants us the luxury of analytical 
treatment and the ability of exploring important characteristics of super-Hubble perturbations and their effect on the 
CMB temperature anisotropy spectrum. As we will see later on in this section, the temperature perturbation arizing 
from a super-Hubble metric perturbation is roughly proportional to the gain or loss in gravitational potential between 
the time of emission of a CMB photon and the time of its observation. This implies that, even if in a more realistic 
case, the super-Hubble perturbation to the gravitational potential decays slightly at late times, this decay is uniform 
at super-Hubble scales. Therefore, we would expect the directional distribution of the temperature perturbation to 
remain the same, modulo a time-dependent factor. Disregarding this time-dependent factor, which is what we are 
doing, makes parameter estimation imprecise, but does not change the qualitative signature of the super-Hubble 
perturbation in the CMB spectrum. 

Given this approximation, the equations of motion for a massless particle now take the form 

Ar) = /i — const (12) 



dA 



z 



A(z)—vt j = (13) 
A- 1 x = v x ,A- 1 y = v v (14) 

The four-velocity = x 1 * is a null vector, which implies 

A- 2 z 2 + v x 2 + v y 2 = if (15) 

Let us define the constant Ao = A(0) = 1 — 2fczo and fix the integration constant ix by setting ix — Ai[2~4|. 
Furthermore, let us fix the affine parameter A by demanding z(0) = 0. This gives 77(0) = 1. We notice that impact 
angles and photon velocity are conformally invariant quantities, and we are therefore free to compute them in any 
conformal frame. Therefore, the physical velocity of the photon in the frame of the observer is u % = A^ 1 x l /r) = i l /Ao, 
and eq. f) 15[) simply states that |u| = 1. We therefore find that the integration constants v x and v y are the velocities 
of the photon in the frame of the observer at the instance of observation. Defining 9 as the observed angle of the 
photon, as measured from the direction of positive 2, we have: 



v T 2 = sin 2 0, cos0 = — u z = z(0)/A o 



(16) 
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The boundary conditions are now defined that allow us to find unique solutions to the equations of motion. Treating 
the case of vertical motion with 8 = or 8 = tt separately, we obtain the following solutions to the equation of motion, 
expressed in terms of the observation angle 9: 

z{X) -\-Xcos8 9 = 0,n (17) 

( r) (9)- ^m(taii(|(0-2fcsin0A))) < 8 < tt 
VW = { L^lnfl- 2fc ™° \) 8 = Q,tt (18) 



2fc cos 6 \ A 



where 770 (6*) = jfc ln(tan(6 l /2)). For simplicity, we have set 77 = at the time of observation. We can now easily 
express the solution in terms of the conformal time r\. We get 




Z<8 Tl) = I 2k V sinS cosh(2fc(r )0 -r,)) 1 ) < 6 < 71 Qg-j 

- ' ^ j^^kn cose _^ e = Qj7r 

Eq. )|19p is parametrized by quantities that are identitical in the two conformal frames, so it is a parametrization 
that also is valid in the physical frame. It allows us to compute the z coordinate of a photon observed at an angle 8 
at any time 77 in its history. Next, we will use this to compute the effect of a super-Hubble metric perturbation on 
the temperature distribution of the CMB. 



B. CMB temperature distribution as measured by a static observer in the CMB rest frame 

In the preceding section, we solved the null-geodesic equations in the background of a super-Hubble metric pertu- 
bation, asserting, for simplicity, that the pertubation is static in comoving coordinates: 

ds 2 = a 2 (77) (-^(z)d?7 2 + A(zy x (dx 2 + dy 2 + dz 2 )) (20) 

where A(z) = 1 + 2k(z — Zq). In this metric, the null geodesic equation for the time coordinate 77 takes the form 
Po = 0, where po = dL/df] = — a 2 At). This implies that p Q is a constant of motion, which we can relate to the energy 

of the photon. The photon energy as observed by a static observer with 4-velocity u^ 1 = ((—.goo) , 0, 0, 0) is 



E = -p^ = -po(-.9oo) 1/2 



Po 



For a photon with observed frequency is, E = hv. We therefore have 



a(j])y/ A(z)v = —po/h = const, 

a generalization of both the gravitational redshift formula of a static universe and the cosmological redshift formula 
of a homogeneous, expanding universe. We will assume that CMB photons were emitted at time 77* at a uniform 
temperature T* . Since the peak frequency of a black-body spectrum is proportional to the temperature of the radiation, 
we can use this relationship to relate the temperature at the time of emission to the observed temperature T{&) of 
CMB photons being observed at angle 8 by a static observer at the origin: 



T(6)y/Ao = T,a(v*WMz(0,V*)), (21) 

where A(z(9,r]*)) defines the gravitational field at the time and place of emission, and z(8,rj^) is given by eq. 
(|19p. Let To be the average CMB temperature at present, and let 0(8) to be the temperature pertubation defined by 
T(8)) = To (l + O (9)). Since To/T* = a(?7*), the redshift relationship gives the following expression for the temperature 
perturbation 0(8) in terms of the gravitational potential at the time of emission for a photon observed at angle 8: 



ew _ (22) 
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With A(z) = 1 + 2k(z — zq) and z(6,rf) given by eq. (|Ti?)) . the expression for the temperature perturbation 9(6*) 
simplifies to 



H "" = ^+1 I (x» i) " 1; (23) 
x * + 1 t/l+ ) X %T^ cos6> 



(x« 2 +i) 

where the constant x* m defined by 2krj* = ln%*. fc is positive while rj* < 0. |?7*| is of the order of the Hq~ , 
the inverse of the current Hubble parameter value, k is the spatial derivative of a super-Hubble perturbation to the 
potential $. Therefore, using the definition of a super-Hubble perturbation, we must have |fcf?*| << 1. Hence, x* is 
slightly less than 1. It is natural then to introduce a small, positive, dimensionless parameter $„ = (1 — x*)/2 that 
will be useful for perturbative expansion of the temperature 0(9). Expanding in terms of (krj*), we have, to lowest 
order in (krj*): 

= - (to/.) + 0((fer/*) 2 ), 

can be interpreted as the loss in gravitational potential since the time of emission for photons observed at angle 
9 = 0. We will refer to <£>* as the potential anisotropy. Then, to second order in $», 9(6*) is 

9(61) = $» cos 9 -$, 2 (\- ^1 + |cos0^cos0^ + O^* 3 ) (24) 

The first-order term <£>* cos 9 equals the temperature perturbation expected from photons following the null geodesies 
of a homogeneous FLRW universe in eq. (|21[) . Hence, higher-order terms stem from deviations from geodesies in a 
homogeneous space. 

In this section, we assumed an observer at rest in the CMB rest frame. In reality, the frame of observation is moving 
relative to the CMB rest frame. Our motion relative to the CMB rest frame has been measured to about 370 km/s 
[l3| . In the next section, we will transform the result of eq. (|23[) to a moving frame of observation. 



C. CMB temperature distribution as measured by a moving observer 

Let O and O' be two inertial observers, O at rest relative to the CMB, while O' moves with constant velocity 
v = vv, where v is the unit vector in the direction of movement and v is the velocity relative to the CMB rest frame. 
Given that O observes a black-body distribution of CMB photons at temperature T(p), with p being the direction 
vector of the observed photon, observer O' will observe a Doppler-shifted black-body spectrum with temperature [HJ 



p ' is 

miii 



T>(p>,v)=T(p) ri V~v% ( 25 ) 
(1 + vp' ■ v) 

is the direction vector of the photon as observed by observer O', and p' ■ v is given by the aberration formula 



P ■ v = --— 26) 

(1 — vp - v) 



Next, let us define two constants 7 = y/ xJ+i an< ^ @ = [i+x* a j ' anows us to write the temperature perturbation 

of eq. (l23l) as 0(9) = „ Z — 1 = n 7 „ « « , where p is the direction vector of the photon, and u is the unit vector 
^ ^ — y v ' vl— p cos y ^l+pp-u 1 c 1 

along the direction of the gradient <9.;<i> of the super-Hubble perturbation, which we hereafter will call the anisotropy 
axis. 

Now, let us use eq. ([25]) to transform the temperature distribution of eq. (|23[) to the frame of observation: 

T\p',v) = - 7 =L=T d (p\v), (27) 
VI + pp-u 



Here, we introduced the dipole distribution Td(p' , v) = Tq\/\ — v 2 /(l + vp' ■ v). Let us first compute p ■ u in the 
frame of observation. Expanding the unit vectors p and u in the CMB rest frame in spherical coordinates as p = 
— (sin 9 p cos 4> p , sin 9 p sin <j) p , cos 9 p ) and u = (sin 9 U cos 4>u , sin 9 U sin 4> u , cos 9 U ) , we get p ■ u = — sin 9 p sin 9 U 



cos 



V 



'V ' 



We can then transform to the frame of observation by using the aberration formula of eq. (|26p 
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and the fact that longitudinal coordinates are invariant under a Lorentz boost along v. We get, using 9 = 6' , <fi = </>' , 
<j) u = <j)' u = and £ = cos#„': 

/(^,^e)=p-u=- 7 - -J— — ((1-u 2 ) Vl-e 2 sin^ cos0 + (£-«) (cos 0-v)) (28) 

(1 — v cos w) (1 — t,v ) V / 

(0, 4>) are now angular coordinates in the frame of observation, and £ is the dot product between the anisotropy 
axis and the direction of motion in the frame of observation: £ = u' ■ v. 

If we assume a completely homogeneous temperature distribution in the CMB rest frame, the distribution observed 
in the moving frame of observation will, according to eq. ([25]) . be the Doppler-shifted dipole distribution Td(p' , v). If 
we average this distribution over the sphere, we get the average temperature in the frame of observation: 

I* S I J dn'T d (p',v) = v^^^To = (l - £ + <V)) To 

We notice that the average temperature in the frame of observation is slightly reduced compared to the static frame. 
The effect is small, though, and proceeding, we will use Tq = To, discarding terms higher than first order in velocity. 



IV. THE EFFECT OF A SUPER-HUBBLE PERTURBATION ON THE CMB POWER SPECTRUM 



Having computed the effect of a super-Hubble perturbation to the gravitational potential on the observed CMB 
temperature distribution, the next challenge is to elicit its effect on the CMB temperature anisotropy power spectrum. 
The standard way of computing this power spectrum is to first remove the dipole induced by the observatory's 
peculiar motion through the CMB sky [16|, then decompose the residual temperature distribution into multipoles, 
and finally compute the covariance matrix of the multipole coefficients. With a homogeneous background geometry, 
this covariance matrix becomes diagonal. With an inhomogenous background geometry, however, the covariance 
matrix obtains nodiagonal entries. Therefore, these nodiagonal entries in the temperature anisotropy covariance 
matrix represent the signature of the inhomogeneity of the background geometry. 



A. Factoring out the Doppler dipole 

Before computing the CMB power spectrum, the CMB dipole must be removed from the temperature data [llij . 
Thereafter, the power spectrum may be computed from the residual temperatures. Our moving observer O' asserts a 
homogeneous background, so we will therefore fit the temperature distribution of eq. (|27|) to that of a Doppler-shifted 
uniform temperature distribution. The procedure is lengthy, but straight forward. We will therefore describe it briefly. 
Let us vary both the velocity and the average temperature by writing the fitted velocity as v" = v + Sv and the fitted 
temperature average as Tq' = Tb(l + ATo). Defining the temperature residual as ST = Td(ft', v") — T'(ft', v), we find 
the best fit dipole by minimizing the averaged squared error J d£!'<5T 2 with respect to variations in velocity v" and 
temperature T^'. The values of Sv and ATo that provide the best fit dipole distribution are 

<^ = $*(l + $*)£-u$*(l-£ 2 ),AT = -i$„ 2 (29) 

Next, we will compute the temperature residual that remains after the Doppler dipole has been removed from 
the temperature distribution. The observed CMB temperature distribution of an assumed homogeneous background 
metric has two factors; the best fit dipole distribution and the observed residual distribution: 

T"(p',v") = T d (p',v")(l + e"(p')) 

Here, 0" is the residual temperature perturbation that remains after factoring out the best fit dipole. 

Similarly, the expected temperature distribution of the inhomogeneous background of eq. (l20l) has three factors; 
the anisotropy factor given by eq. (|23p , the dipole distribution and some arbitrary residual temperature distribution 
T = I + given a priori: 

T\f>', v) = 1 up', v) (i + e(p')) 

y/l + Pf{p',0 

Using eq. (|2~T|) , we may then relate the observed temperature pertubation 0" to the given perturbation 0: 

Q"(p') = —= 7 V }. (1 + Q(p')) - 1 (30) 

Eq. (|30p will be used later in the final computation of the CMB anisotropy power spectrum. 



B. CMB anisotropy spectrum 

The temperature distribution Q"(j>) can be expanded in terms of spherical harmonics Y l m (p) and multipole coeffi- 

I 

cients a\ m as 0"(p) = S £ ai m ^ m (p)- This equation can be inverted, expressing the multipole coefficients a\ m in 

/ m=—l 

terms of the field Q": 

a lm = J dnY l m (p)*e"(p) (31) 

Let us define the CMB covariance function as the expectation value of the square of the multipole coefficients: 

C z"fc" = ( a fc« a L) (32) 

For a homogeneous background geometry, C^jf 1 can be expressed in terms of the Ci of the conventional CMB 
anisotropy power spectrum, in which case it takes the simple, diagonal form 

Now, let us compute C;™'™ for the inhomogeneous background of eq. . Using eq. (|31[) , we get 



c# n = Jdti'J <mY k n (p>)Yr(py < e"(p')e"(p) > 

Now, Q" is given by eq. ([3H| . Define A0 as 

7 Td , (f3>) , ^ (33) 

Now, using eq. (|30j) , the observed temperature perturbation 0" can be expressed in terms of A0 as follows: 

e'V) = (i + A6) (l + e(p')) - 1 

This allows us to express Cj™' n in terms of A0: 

c$ n = | d«'r fe "(p')Ae(p')* | dnYr(pfAe(p)+ J dn'Y^ip') (i + A©(p')) / dm^)* (i + AG(p)) < ©(p')*e(p) > 

(34) 

The expectation value < Q(p') &{p) > can be expressed in terms of the Fourier-transformed temperature distribu- 
tion 0(fe, k ■ p), the matter overdensity S(k) and the matter power spectrum P{k) (see eq. (8.65) of ): 

d 3 k @(k,k-p')* e(k,k-p) 



<e(p')*e(p) >= f -^p(fe) 

J (Z7r) 



(2tt) j v ' S(k) S(k) 

The first term on the right-hand side of eq. (|34[) is second order in the perturbation variables and v, so we will 
discard it in the following. This gives 

d 3 k &(k,k-p')* Q(k,k-p) 



c% n = J dn% n (p>)(i + AG(p>)) J dnYr(pT(i + Ae(p)) J ^4^*) 



(2tt) 3 v ' S(k) 5(k) 



First, we may expand 0(fc, k ■ p) in terms of multipoles 0;(fc): 



0(fc, k-p) = {21 + 1) Vi{k ■ p)e l (k) 

Furthermore, the Legendre polynomials Vi can be expanded in terms of spherical harmonics: 

ZL + 1 m=—l 



Now, define the parameter r = — £ 2 . The covariance function C;™'™ then takes the form 

C^=<34,*<W + T(C fc + Q) <r (35) 
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where C\ is defined by 

Q = - f dkk 2 ^\e t (k)\ 2 (36) 
7T J 5(k) 

and the coefficients K7Y 1 are defined by the integral 

K%» = r- 1 J dnAe(p)Yr(pTY h n (p) (37) 

We see that Ci represents the variance of the multipole coefficients a\ m , and its definition in eq. ([55)1 is indeed the 
conventional formula for computing the variance of the multipole coefficients a\ m for a homogeneous background, see 
p. 242 of [12| . Notice that when the expression for C ; ™'™ in eq. (|33|) was computed, terms of order r 2 were discarded. 

Eq. (|35p is the general form of the covariance function C ; ™'™. If the coefficients K™^™ are nonzero, the covariance 
function is nondiagonal. However, we have yet to compute the coefficients K, £ and prove that there are indeed 
nondiagonal entries in the covariance function. 



C. Nondiagonal entries of the covariance function 

Next, we will compute the nondiagonal coefficients Kj n ^ n of the covariance function C ; ™'™. The coefficients are given 
by eq. (|3"Tj) . Before we are able to compute the integral, we must compute A0 from its definition in eq. (f3"3")) . To first 
order in the perturbation variables and v, we get 

A 6 = $* y/l -£ 2 sin 9 cos 4> (38) 
Spherical harmonics may be written in terms of the associated Legendre functions P, m (see p. 670 of [17]): 



Y r{0A) = ^ 2l l 1) l + ^ ] pr{^ey^ 

■Kj k' n can then be expressed as follows, defining the angular variable /i = cosQ: 

/ (21 + 1) (Z - m)7 / (2k + 1) (fc ^n)\ ) r J i{n _ m)(f> 

V ' -1 o 

Now, we may compute the integral over the longitudinal variable <p: 

/d^-^cos^ = ir(S m , n+1 + 5 m+hn ) 
o 

Furthermore, we may use the following recursion relationship for the associated Legendre functions fl8l Il9j: 

y/T^PT- 1 ^) = 727TT) - p ™i(^)) 

as well as the orthogonality relationship (see p. 66 of [19] ): 

-l 

By introducing the coefficients 



kkm. —Km<l 

2l+1 S m S (3Q) 
Iml > / 
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we obtain the following expression for the coefficients K™^ n : 

Kk" = \ ((cTcTSrn.n+l - df™ C^™ 5 m , n - X ) 5l,k+l + (c?c£<5 m+ l,„ - CpCfc "<5 m -l,„) (40) 

The coefficients K7u n have the following symmetries: 

«IS n = *IT ( 41 ) 

and 

^ n = -^T ,_W ( 42 ) 

D. Estimating the potential anisotropy vector 

The perturbation parameter <£>* measures the loss in gravitational potential for photons arriving along the anisotropy 
axis v! . Therefore, measures the gain or loss in gravitational potential for photons arriving along the dipole axis 
at = 0, while $»a/1 - £ 2 measures the gain or loss in gravitational potential for photons arriving from directions 
orthogonal to the dipole axis, at 9 — tt/2. 

We refer to the vector as the potential anisotropy vector. yl — £ 2 is then the transverse component of the 
potential anisotropy vector, which is the component orthogonal to the dipole axis, while is the longitudianal 
component of the potential anisotropy vector, which is the component along the dipole axis. Define A = The 
parameter t = $*yl— £^ has already been introduced. We will refer to r as the transverse potential anisotropy 
and A as the longitudinal potential anisotropy. In order to estimate the direction of the the anisotropy vector, we 
also need to estimate the longitude of the anisotropy axis in galactic coordinates. Alternatively, by noticing that the 
temperature anisotropy A0 = r sin9 cos</> has extrema at 8 = 7r/2,</> = ±7r, we could estimate the direction of the 
potential anisostropy based on an estimate of the direction of maximal temperature anisotropy perpendicular to the 
CMB dipole axis plus estimates of the parameters r and A. If w is the observed direction of maximal temperature 
asymmetry, we have 

= tw + Xv (43) 

We see from eq. (|40[) that to lowest order in the perturbation parameters, we are only able to estimate the transverse 
potential anisotropy r. In order to estimate the longitudinal potential anisotropy, the covariance function C ; ™'™ must 
be evolved to higher order in the perturbation variables. In the scope of this paper, we will restrict ourselves to lowest 
order estimates, so we will therefore only be able to estimate the transverse potential anisotropy. 

From eq. (|4"0)) we obtain four sets of estimators for r, one for each of the four terms. The first and third terms give 
the same estimator, and so do the second and fourth terms. The remaining two estimators can be identified by using 
the symmetry relation of eq. (I42|) . The basic estimators can then be written 



/L "M-1 



(I -2),.. .,(1-1) (44) 



This basic estimator can be used for estimation of the transverse potential anisotropy r from multipole coefficients 
ai m derived from CMB observations. For each multipole index I, we may sample r 2(1 — 1) times. If we choose a 

N+l 

set of N multipoles ranging from I = 2 to I = N + 1, we will therefore be able to retrieve S 2(7 — 1) = N(N + 1) 
independent samples of r, thus reducing the uncertainty of the estimate considerably. 



E. Hemispherical power asymmetry 

Finally, let us have a glance at the hemispherical power asymmetry that would result from the temperature 
anisotropy A9 = rsin^cosi/) of eq. (|3"8|) . Designating (9 = ~,(f> = tt) as the north pole, the northern hemi- 
sphere is parametrized by 7r/2 < <j) < 37r/2. We would therefore expect a maximal temperature difference between 
the two hemispheres to be ~ 2t0, yielding a power difference between the two hemispheres of the order of rC;. 
Furthermore, given that the north pole is at 9 = 7r/2, the axis of maximal power asymmetry should be perpendicular 
to the CMB dipole axis. Given a CMB dipole at (l,b) = (264°, 48°) [13j and a best fit power asymmetry axis at 
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TABLE I: Separation angle between the dipole axis and the axis of maximal power asymmetry for different results of Eriksen 

et. al. [6J. a) small sky cut applied, b) large sky cut applied 

Data (l,b) Separation angle 

ILC a (225°, -27°) 83° 

ILC 6 (208°, -27°) 90° 

Q-band (222°, -35°) 91° 

V-band (205°, -19°) 85° 

W-band (204°, -31°) 96° 



(I, b) = (225°, — 27°)[a|, the angle separating the two axes is 83°. Table U shows the angle of separation for each of 
the 5 different analyses made by Eriksen et. al @. In each case, the angle of separation between the CMB dipole axis 
and the axis of maximal power asymmetry lies between 83° and 96°, consistent with the prediction by our model. 

We notice that eq. ([4"3")) can be applied to estimate the direction of the potential anisotropy from estimates of the 
power asymmetry axis w, given estimates of parameters r and A. 

An interesting side note is that the anomalous CMB Cold Spot (20l - [22l | is located in the same region of the sky as 
the north pole of the power anistotropy axis: The cold spot is located at (I, b) = (209°, —57°). Given that there is a 
CMB power asymmetry with a minimal power in this particular region and assuming this power asymmetry can be 
explained by the presence of a super-Hubble inhomogeneity in the gravitational field, the cold spot will be less extreme, 
because the expectation value for the CMB temperature would be lower in this particular region. Consequently, the 
cold spot should be less anomalous in this CclSG, els already noted by Eriksen et. al. [6||. 



V. CONCLUSIONS 

Prompted by the puzzling evidence of an hemispherical power asymmetry in the CMB Q, this paper started with 
two simple questions: 1): How could a large-scale inhomogeneity in the gravitational field of super-Hubble size - 
larger than the observable universe - be observed? And 2): Could an hemispherical power asymmetry in the CMB 
be caused by a super-Hubble scale inhomogeneity in the gravitational field enclosing the present Hubble volume? 
Posing these questions, we made no assumptions about the origin of this perturbation or any physical mechanism 
that might have caused it. Our goal was simply to determine what effect the presence of such a perturbation to the 
gravitational potential would have on the observed CMB temperature and to elicit its observational signature in the 
CMB spectrum. 

In order to tackle these questions by analytical means, we introduced an idealized model in which large-scale 
perturbations to the gravitational potential stay constant in comoving coordinates. This model is an unrealistic ap- 
proximation at all but the largest scales. At large, super-Hubble scales, still being a very crude approximation[25j, 
it is nevertheless a useful idealization that grants us the luxury of analytical treatment and the ability of explor- 
ing important characteristics of super-Hubble perturbations and their effect on the CMB temperature anisotropy 
spectrum. We find that the temperature perturbation arizing from a super-Hubble metric perturbation is roughly 
proportional to the loss in gravitational potential between the time of emission of a CMB photon and the time of its 
observation. Therefore, even if in a more realistic model, the super-Hubble perturbation to the gravitational potential 
decays slightly at late times, this decay is uniform at super-Hubble scales. Therefore, given a more realistic model, 
we still expect the directional distribution of the temperature perturbation to remain the same as in our idealized 
model, modulo a time-dependent factor. Disregarding this time-dependent factor, which is what we are doing, makes 
parameter estimation imprecise, but does not change the qualitative signature of the super-Hubble perturbation in 
the CMB spectrum. 

Using our idealized model, we were able to solve the null-geodesic equations of motion and determine the resulting 
CMB temperature distribution in the CMB rest frame. Transforming the resulting temperature distribution to a 
moving frame of observation and removing the best fit Doppler dipole, we obtained a residual temperature distribution. 
By expanding the residual temperature distribution into multipoles, we obtained our main result: the covariance 
function of the multipole coefficients. 

The covariance function contains, in addition to the C; entries of the conventional CMB temperature anisostropy 
power spectrum, non-diagonal entries. We find that the Ci entries of the conventional anisotropy power spectrum are 
insensitive to the strength of the potential anisotropy. Thus, the non-diagonal entries, which are not present in the 
case of a homogeneous background geometry, constitute the main signature of a large-scale inhomogeneity in the 
background geometry of the universe. This anwers the first of our two initial questions. 
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Regarding the second of our initial questions, an inhomogeneity in the gravitational potential of super-Hubble size 
would yield a power asymmetry in the CMB with maximal asymmetry at an angle of 90° to the CMB dipole axis. 
The power asymmetry that was observed appears to be at an angle that lies between 83° and 96° to the CMB dipole 
axis, which is consistent with the prediction of our model. This is suggestive of a simple explanation for the CMB 
power asymmetry, because it implies that the location of the observed power asymmetry in the CMB sky can be 
accounted for by a large-scale inhomogeneity in the gravitational field enclosing the present Hubble volume. At this 
point, we will not claim that our model can completely account for the CMB power asymmetry. More work remains. 
In particular, it still remains to be seen whether this model is also able to account for the power asymmetry in 
quantitative terms. 



Appendix A: Super-Hubble solution to the first order perturbation equations for ACDM 

Let us start by taking a look at the evolution equation for the temperature perturbation O = AT/T. We use 
conformal Newtonian gauge and apply the notation and exposition of Dodelson (see ch. 7 of [l2|). is a function, 
not only of conformal time rj and position x, but also of photon direction, p l . The first order perturbation equation 
for 0(?7, x,p) is 

dQ d$ -> 

■7T- + 1T +fd i e + p l d i ^ = n e a T a(& -&+p- v b ) (Al) 
or] orj 

where n e is the number density of free electrons, vi, is the baryon velocity and ot is the Thomson cross section 
for electron-photon scattering. 0o is the temperature monopole: ®o(r),x) = J dil<d(rj,x,p). Seeking super-Hubble 
solutions to eq. (|Aip . we may discard the terms containing spatial derivatives. Integrating the equation by j- J dQ, 
eq. (I All) simplifies, leaving an equation for the monopole O and the potential <I>: 

— —-ii (A2) 
drj drj 

Here, we have used that f dQp ■ u = for any direction-independent vector u, which implies that J dVtp ■ ^u^ = 0. 

For super-Hubble perturbations we may therefore discard baryon and photon interactions, which occur at much 
smaller scales. At super-Hubble scales, perturbations to the dark matter and baryon distributions can be treated as a 
common distribution of matter. Similarly, at super-Hubble scales, perturbations to the distributions of photons and 
neutrinos can be treated the same, as a common distribution of relativistic radiation. 

The two Einstein equations determining the evolution of the metric perturbations take the form 

m^- - d 2 $ - m 2 ^ = iirGa 2 (p$6 + 4 P ^e ) (A3) 
arj 

d 2 ($ + *) = (A4) 

Here, S = 8p m /pm is the matter overdensity, while 0o is the monopole moment of the radiation distribution. p$ 
is the zero-order matter density, which takes the form p$ = Pcr^m/i 3 , where p CI is the critical density today and 
£l m is the present matter ratio. pi°^ is the zero-order radiation density, which takes the form p 7 . = p cl -Q r /a' 1 . 

The equation for the dark matter density perturbation is 

85 5$ 

_ +3 _ + ai „- = 0, (A5) 

while the equation for the dark matter velocity v l is 

f) i 

-H-+Hv i + d i V = Q (A6) 

OTj 

Eq. (|A4[) is easily solved by setting ^ = — $. Again, as we seek super-Hubble solutions to the equations, we discard 
terms with spatial derivatives. In that case, the velocity equations for matter decouple from the other equations, and 
we are left with the three equations for the three perturbation variables $, 0o and 6. Furthermore, as the equations 
now only contain time derivatives of these variables, we may separate each perturbation variable X(rj,x) into a 
time-independent factor and a time-dependent factor. The time-independent factors can be determined by applying 



13 



initial conditions. We apply adiabatic initial conditions, which are <f>(i]i,x) — $>q(x), ®o(r]i,x) = ^$o(x) and 
S(rji,x) = ^<f>o(x), where r\i is an arbitrarily chosen initial time. &o(x) is an arbitrary perturbation in the potential 
satisfying |9i$o| << H and 9 2 $o << T~L 2 - Hence, the time-indepdent factor of each of the three perturbation 
variables are &o(%)- We will therefore rewrite the variables as follows: <do(r],x) — > Qq(t])$>o(x), 5(r],x) — > S(rf)$>o(x) 
and $(ry, x) — > §>(r))$>o(%)- This allows us to factor out the spatial dependence of the equations entirely, and we are 
left with the following ordinary differential equations, using overdots to denote differentiation with respect to r\: 

6 + $ = (A7) 
3H$ + 3-H 2 $ = 47rGa 2 (p(° ) <5 + 4^9 ) (A8) 
5 + 3$ = 0, (A9) 

Applying adiabatic initial conditions = 2Qo(r]i) and S(rji) = 3Qo(r]i) at an initial time rji, allows us to integrate 

eqs. (J7| and (|A"9|) . We get 

Mv) = 8 -^-Ari) = \-mri) 

Changing to a new temporal variable x = In a, eq. (|A8|) can now be written 

d£ _ 1 {n m e' 3x + ln r e- 4x ) /9 
~dx + 2 £l\ + n m e~ 3x + n r e- 4x V 2 



' on ,0 „-3z 1 o B -4x I n I \-"-^ u J 



There are no known analytical solutions to eq. (|A10[) , but we may find approximate solutions that cover the entire 
history of the universe. First, let us seek an approximate solution that covers the radiation and matter-dominated 
eras. In this case a 3 << 1. Introducing a new variable y — jf^e x — jf-a allows us to rewrite eq. (|A10[) as follows: 

U dy \ 2y + l) 4(y + l) V ' 

Eq. (|A11|) has the following solution that satisfies the initial condition that requires $(y) — > 1 as y — » 0: 

^ = (-16- 8y + 2y 2 + 9y 3 + 16^/1+^) (A12) 

We see that $ has an almost constant value of 1 throughout the radiation era, then drops slightly to approach 9/10 
in the matter-dominated era. 

Next, let us seek an approximate solution to eq. (|A11[) that is valid in the matter-dominated era and later. In eq. 
(jAlip . we can now discard the radition terms. Defining a new variable z = j^e 3x , we get the equation 

d$ 2z + 5 3 1 

Z— + —, -T$ = T " 



dz 6 (z + 1) 4^ + 1 
Its solution can be expressed in terms of the hypergeometric function 2-F1: 



* = (1 + ^)2^(1,1;^;-^) (A13) 



3' 6 

This solution has the value $(0) = 9/10. At present times, <I> ~ 0.7 for Q m = 0.3. 
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